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ABSTRACT 

Context. Observations of the atomic and molecular line emission associated with jets and outflows emitted by young stellar objects 
provide sensitive diagnostics of the excitation conditions, and can be used to trace the various evolutionary stages they pass through 
as they evolve to become main sequence stars. 

Aims. To understand the relevance of atomic and molecular cooling in shocks, and how accretion and ejection efficiency evolves with 
the evolutionary state of the sources, we will study the far-infrared counterparts of bright optical jets associated with Class I and II 
sources in Taurus (T Tau, DG Tau A, DG Tau B, FS Tau A+B, and RW Aur). 

Methods. We have analysed HerschelfPACS observations of a number of atomic ([Oi]63/mi, 145pm, [Cn]158/jm) and molecular 
(high-J CO, H 2 0, OH) lines, collected within the Open Time Key project GASPS (PI: W. R. F. Dent). To constrain the origin of the 
detected lines we have compared the obtained FIR emission maps with the emission from optical-jets and millimetre-outflows, and 
the measured line fluxes and ratios with predictions from shock and disk models. 

Results. All of the targets are associated with extended emission in the atomic lines; in particular, the strong [O i] 63 fim emission 
is correlated with the direction of the optical jet/mm-outflow. The line ratios suggest that the atomic lines can be excited in fast 
dissociative J-shocks occurring along the jet. The molecular emission, on the contrary, originates from a compact region, that is 
spatially and spectrally unresolved, and lines from highly excited levels are detected (e.g., the o-FLO 8 ig - 7 7 line, and the CO J=36- 
35 line). Disk models are unable to explain the brightness of the observed lines (CO and FLO line fluxes up to 10 _15 -6 10 -16 WrrT 2 ). 
Slow C- or J- shocks with high pre-shock densities reproduce the observed H 2 and high-J CO lines; however, the disk and/or UV- 
heated outflow cavities may contribute to the observed emission. 

Conclusions. Similarly to Class sources, the FIR emission associated with Class I and II jet-sources is likely to be shock-excited. 
While the cooling is dominated by CO and ILO lines in Class sources, [O i] becomes an important coolant as the source evolves and 
the environment is cleared. The cooling and mass loss rates estimated for Class II and I sources are one to four orders of magnitude 
lower than for Class sources. This provides strong evidence to indicate that the outflow activity decreases as the source evolves. 

Key words. Astrochemistry - Stars: formation - ISM: jets and outflows - ISM: molecules - ISM: individual objects: T Tau, DG Tau 
A, DG Tau B, FS Tau A, FS Tau B, RW Aur 



1. Introduction 

Theor etical models (e.g., IShu et all Il994t iKonigl & Pudritzl 
2000) predict a tight correlation between the accretion of 
matter onto a young star and the ejection in winds and/or 
jets. Measurements of stellar accretion and mass loss (e.g., 
lHartigan et alj [l995l) support the general picture presented by 
these models, but the uncertainties of these measurements are 
too large to provide a quantitative test of the predictions of the 
ratio of the mass accretion rate to the mass loss rate. Sources in 
the earliest stages in their evolution (Class 0) are not visible, and 



* Herschel is an ESA space observatory with science instruments 
provided by European-led Principal Investigator consortia and with im- 
portant participation from NASA. 



are often indirectly identified by means of their strong ejection 
activity, which is manifested in the form of bipolar parsec-scale 
molec ular outflows of ten observed at millimetre wavelengths 
(e.g., Bachillerl ll996l) . The ejection associated with evolved, 
optically visible T Tauri stars (i.e., Class II) is instead usually 
traced by bright blue- and red- shifted forbidden emission 
lines p resent at optical and near-infrared (NIR) wavelengths 
(e.g., lHartiganetai1ll995l) . For Class I sources and some 
Class II sources both th e molecular outflow and t he optical 
jet ha ve been observed (Gueth & Guilloteau 1999; iPetv et all 
2006). These observations show that the two components are 
connected. The optical/NIR forbidden lines trace hot (~10 4 K) 
atomic gas, which is believed to have been extracted from the 
disk, and accelerated in the observed fast and collimated jets 
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(velocities up to hundreds of kms~' and jet widths smaller 
than 200 AU). The millimetre observations, instead, trace cold 
(~10-100 K) and slow (tens of kms -1 ) gas, which is thought to 
be ambient gas that has been set into motion by the jet propa- 
gation (i.e., jet-driven m olecular outflow, e.g., Raga & Cabrill 
119931; ICabrit et all Il997t) . However, colli mated high velocity 
molecular gas (velocity up to ~60 kms -1 , Lefl och et al.l 12007) 
has also previously been detected at millimetre wavelengths, 
questioning this simple picture, and suggesting that molecules 
can also be extracted from the disk and acce lerated in the jet 
dPontoppidan et alj|201 UlPanoglou et"aTll2012l) . 

In this context, observations at far-infrared wavelengths al- 
low us to trace the intermediate warm gas component in the 
jet/outflow system. Prev ious observations from the Infrared 
Space Observatory (ISO. iKessler et all 1 19961) targeting outflow 
sources have shown emission in a large number of atomic ([O i], 
[Cn]) and molecular (H2O, CO, OH) lines. Despite the very 
low spectral and spatial resolution offered by ISO, analysis of 
the line fluxes and ratios indicates that the bulk of the detected 
[Oi] and molecular emission is most l ikely to be excited in the 
shocks occurring alo ng the jet/outflow (Nis ini et alJ [T996. 1999; 
iGiannini et al.ll200~lb . The line fluxes have previously been used 
to estimate the cooling in the atomic and molecular lines, and 
to quantify the outflow efficiency as the ratio between the to- 
tal luminosity radiated away in the far-in frared lines, L (FIR), 
and the source bol ometric luminosity, Lj, i (IGiannini et alfc OOl ; 
iNisini et"a"Dl2002l) . However, because of its limited sensitivity, 
ISO observations have been restricted to studies of bright and 
extended outflows from Class and I objects. The ESA Herschel 
Space Observatory (Herschel) has allowed, for the first time, ob- 
servations of the FIR counterparts of optical jets associated with 
Class I and Class II sources whose environment has been largely 
cleared. 

As the source evolves, the accretion/ejection activity is ex- 
pected to decrease, with the surrounding cloud material being 
either accreted onto the star, or dispersed by the jet. As a con- 
sequence, the optical jet will become visible while the emission 
at far-infrared wavelengths should be expected to be fainter and 
less extended than in Class sources. However, FIR observa- 
tions of the ejection activity associated with more evolved Class 
I and II sources is interesting for the following reasons: 

- The source, disk, and accretion properties of T Tauri stars are 
well-known, and so we can study the correlation between the 
detected ejection phenomena and these properties. Thus, we 
can estimate the mass ejection to mass accretion ratio; 

- It is usually the case that Class sources are observable only 
at millimetre wavelengths, whilst the ejection activity from 
T Tauri stars is detected only in the optical. FIR lines can 
however be detected from Class to Class II sources, there- 
fore facilitating a way to form an evolutionary picture of jet 
activity; 

- The detection of molecular emission in sources whose en- 
vironment has been cleared may support the idea of a disk- 
wind molecular component pro viding strong constrai nts to 
existing models of jet launching (Panoglou et alj|20"l2h . 

Thus we have analysed the FIR emission from five Class I 
and II sources in Taurus (d~140 pc). These sources were ob- 
served as part of the Herschel Open Time Key project GASPS 
(GAS in Protoplanetary Systems , PI: W. R. F. Dent) using the 
PACS inte gral field spectromet er dPoglitsch et al.l2010h on board 
Herschel (Pil bratt et al.l 12010b . These PACS observations pro- 
vide maps of the emission in a number of atomic and molec- 



ular lines (see. iMathews et a"l1l2010l and Dent et al., in prepa- 
ration). The five sources presented in this paper form a subset 
of the Taurus sample analysed in Howard et al., in preparation. 
These sources were selected on the basis of their association 
with bright and extended stellar jets detected in t he typical [SnL 
JN ill , and IQil optica l forbidden lines (e.g., iHartigan et al.l 
11995b iHirth et al.lll997l) . The details of the observations, and the 
applied data reduction processing, are described in Sect. [2] In 
Sect. [3]we show the obtained spectra and maps and compare the 
spatial distribution of the atomic and molecular emission with 
that of the associated optical jets. In Sect.|4]we will compare the 
observed line fluxes and ratios with predictions from both disk 
and shock models. Hence, we will use the observed line fluxes 
and the results from shock modelling to estimate the far-infrared 
cooling and the mass loss rate. The comparison with the values 
estimated for Class and I sources will allow us to place Class II 
sources into an evolutionary picture. Finally, in Sect.[5]we sum- 
marise our conclusions. 

2. Observations 

The observations analysed in this paper have been acquired us- 
ing the Herschel-PACS integral field spectrometer. The criteria 
used to select the sources are explained in Sect. 12. II Information 
on the instrumental settings and on data reduction processing 
are given in Sect. 12.21 and then in Sect. I2.3l we will explain the 
procedure to distinguish extended and unresolved emission as 
observed with the PACS spectrometer. 

2.1. Sample selection 

Many of the sources in Taurus are associated with stellar jets 
based on the detection of extended emission in optical forbid- 
den lines, such as [Sn], [Oi], and [Nn] lines in the 6300-6700 
A wavelength range, which are thought t o be excited in the 
shocks occurrin g along the jet (see, e.g., Hart igan et al.lll99"5t 
IHirth et aT1ll997l) . Within the Taurus sample (Howard et al., in 
preparation), we selected five well-known jet sources (T Tau, 
DG Tau A, DG Tau B, FS Tau A+B, and RW Aur) having the 
following characteristics: 

- they are associated with bright jet emission in optical forbid- 
den lines, extending on spatial scales larger than the PACS 
spatial sampling (i.e., ~9'.'4); 

- they show extended emission in the [Oi] 63 pm line (see 
Sect.OandOfor details); 

- they also show emission in the atomic [Oi] 145 pm and 
[Cn] 158 pm lines and in a number of H2O and high-J CO 
lines (see Sect. which make it possible to compare ob- 
served line ratios and fluxes with predictions from disks and 
shocks models (see Sect.[4]i. 

These characteristics make these sources ideal candidates to 
study the far-infrared atomic and molecular counterpart of op- 
tical stellar jets. Note that, with the exception of RW Aur, their 
environment is not completely cleared: T Tau, DG Tau B, and 
FS Tau B are assoc iated with CO outflows detected at millime- 
tre wavelengths (lEdwards & Snehl 119821; iMitchell etafl 119971; 
iDavis et aLll2010l) . and arc-shaped reflection nebulae associated 
with the optical jet have been observed for T Tau (W-E jet, blue 
lobe), DG Tau A (blue lobe) , and DG Tau B (both lobes) by 
Stapelfel dTet al.l (Il997lll998h . This suggests that these nebulae 
are illuminated outflow cavities. 

With the exception of T Tau, which has previously been ob- 
served in the FIR with ISO showing emission in [0 1], [C n], OH, 
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Table 1. Source properties: stellar and bolometric luminosity, SED class, and position angle of the associated optical jet. 



Source 


L, 




Class 


PA; a 


Ref 




(Lq) 


(L ) 




o 




T Tau (N, Sa+Sb) 


7.3,- 


15.5, 10 


II, I 


180, 270 


(1) 


DG Tau A 


3.2 


6.36 


II 


226 


(2) 


DG Tau B 




-1.1 


I 


115 


(3) 


FS Tau (Aa+Ab, B) 


0.15+0.17, - 


1.4, >0.5 


II, I 


-, 55 


(4) 


RW Aur (A+B) 


1.7+0.4 


3.2 


II 


120 


(5) 



Refer ences. (1 ) L» (only for T Tau N ) by lWhite & Ghezl d2001h. L bo i bylKenvon & Hartmam] dl995T) . Class bvlFurlan et al.l j2006h:lLuhman etaD 
d2010l).PA,, f by lSolf & Bohm| dl999T) : (2) L» and Class bv lLuhman etall(l2010h:lRebull et alJd2010t> . U„, bylKenvon & Hartmanrtd 19951). PA, e , b y 
Mundt & Fried ( 1983); (3) U„i average value of the range estimated bv lKruger et al] {201 1[). Class by Luhman et al. (2010); Rebull et al. (2010), 
PA ; „ by iMundt & Fried {1983); (4) L,, h boh and Cla ss for FS Tau Aa +Ab by lHartigan & Kenvonl d2003l). IKenvon & Hartmanr] d!995l). and 
ILuhman et al.l d2010h. L M and Cl a ss for FS Tau B bvlStark eTall d2006h. andlLuhman et al.l d2010h: IRebull et al.1 d2010h. PA, „, bvlMundt et al] 
d!984l); (5) L, by IWhite & Ghezl d2001l) . h bo , by IKenvon & Hartmard ( fl995l) . Class by iFurlan et al] d2006T) ; ILuhman et al.l d2010l) . PA,„"bv 
lHirthetal](fl994tl 19971) . 



CO, and H2O lines (Spinoglio et al . 2000), and a tentative detec- 
tion o f the [0 1] 63 /mi line at 3.5cr for DG Tau A (ICohen et al.l 
1988), the sources in the sample have not previously been ob- 
served at FIR wavelengths. 

The stellar luminosity, bolometric luminosity, and Class of 
the selected sources are summarised in Table [TJ along with 
the position angle of the associated optical jets. The class has 
been estima t ed by means of recent Spi t zer observations by 
IFurlan et al.l d2006l) : ILuhman et all d201(H) ; IRebull et al.1 d2010b 
from the source spectral energy distri bution (SED) in the in - 
frared accor ding to the classification of lLada & Wilking| (ll984): 
lLadal dl987l) . For multiple systems (T Tau, FS Tau, RW Aur) we 
list the stellar, bolometric luminosity, and the class of each re- 
solved component, along with the PA of the jet associated with 
each of them. However, these systems, and their associated mul- 
tiple jets are not resolved with PACS. Thus the continuum and 
line flux values quoted later in this paper refer to the whole sys- 
tem, as explained in the following section. 



2.2. Instrumental setting and data reduction 

The observations were acquired between February 2010 and 
March 2011. The PACS integral field unit (IFU) allows simul- 
taneous imaging of a 47"x47" field of view, resolved into 5x5 
spatial pixels of 9.4" x9.4" (also called spaxels). In spectro- 
scopic mode for each spaxel a ID spectrum is recovered simul- 
taneously in a selected spectral range in the blue (B: 51-105 /mi) 
and in the red (R: 102-220 //m) arm of the spectrometer. The ob- 
servations were carried out in the chop-nod mode to remove the 
background emission and with a single pointing on the source. 
Spectroscopic observations in line mode (PacsLineSpec) were 
obtained by performing one nod cycle for a total on-source inte- 
gration time of 1 152 s. These simultaneously covered a selected 
wavelength range in the blue and in the red arms. Spectroscopic 
observations in range mode (PacsRangeSpec) were acquired by 
performing 1 nod cycle for T Tau, DG Tau A, and RW Aur, and 
2 nod cycles for DG Tau B, and FS Tau A+B. These scanned 
three wavelength ranges in the blue, and simultaneously, in the 
red arm, for a total integration time per spectral range of 1592 
s (1 nod cycle) and 3184 s (2 nod cycles). In some of the ob- 
servations: specifically in the case of DG Tau and DG Tau B, the 
target was not centred on the central IFU spaxel, but at a position 
~6.7" away from it. The order, arm, spectral coverage, spectral 
resolution, integration time, and targeted lines relative to the ac- 
quired line and range spectra are summarised in Table [2] The 



observation identifiers (OBSIDs) are summarised in Table |A~T1 
in Appendix [A] 

All the data were reduced using HIPE 4.0.1467. The PACS 
pipeline included corrections for: saturated and bad pixel re- 
moval, chop subtraction, relative spectral response function cor- 
rection, flat fielding, and mean of the two nods. The spectra were 
Nyquist binned in wavelength, with non-overlapping bins set to 
half the width of the instrumental resolution. An aperture cor- 
rection was not applied because (i) in some cases the target was 
not centred onto the central IFU spaxel and (ii) extended emis- 
sion was detected in some of the detected lines. The continuum 
flux was recovered by summing the emission over the 5x5 ar- 
ray, then by applying a first-order polynomial fit to estimate the 
continuum level at the line rest wavelength. The error on the con- 
tinuum was calculated as the standard deviation of the difference 
between the estimated continuum level and the observed one in 
a region from 2 to 10 instrumental FWHM from the line rest 
wavelength. The line flux was recovered by summing the emis- 
sion over those spaxels where the line wass detected to avoid the 
suppression of faint lines by the noise in the outer spaxels. The 
integrated flux of the detected lines was estimated by measur- 
ing the total flux under the gaussian line fit. We calculate the lcr 
flux by integrating a gaussian with height equal to the continuum 
RMS, and width equal to the instrumental FWHM; we report 3<x 
upper limits for non-detections. All of the spectral ranges cov- 
ered by our observations in the spaxel where the line emission is 
maximum are shown in Fig. [TJ and the continuum and integrated 
line fluxes are summarised in Tables[3]and|4]. 



2.3. Extended and unresolved emission 

Fig. |2](left panels) shows the ID spectra in the 63. 12 - 63.25 /mi 
range for the 25 spaxels of the PACS array of the 5 sources in 
our sample. All of the maps show bright [0 1] 63 /mi emission 
in a number of spaxels. This suggests that the line may origi- 
nate from an extended region such as the jets detected at opti- 
cal wavelengths, rather than from a circumstellar disk. Indeed, 
if the spectroscopic point-spread-function (PSF) is smaller than 
the spaxel size, the emission from the source and the circumstel- 
lar disk (typical disk sizes in Taurus are <500 AU) will be spa- 
tially and spectrally unresolved with PACS, and thus confined 
in the central spaxel of the 5x5 integral-field array. On the con- 
trary, if the emitting region is more extended than the spaxel size 
(9'.'4x9'.'4) the line and continuum may be detected in the outer 
spaxels. 
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Table 2. Spectroscopic Herschel/PACS observations: order, arm, spectral coverage, resolution, integration times, and targeted lines. 



Order 


Arm 


Spectral Coverage 


K 


1 int 


Lines 






Qan) 


(kms ') 


(s) 




i 


D 


OZ.yj - 03.43 


on 
65 


1 1 CO 

1 Ijl 


|OiJ Pi- P 2 , o-H 2 818-/07 


1 


R 


180.76 190.29 


200 


1152 


DCO + J=22-21 


2 


B 


71.82-73.33 


162 


1592, 3184 


o-H 2 7 07 -6 16 , CH + J=5-4, CO J=36-35 


1 


R 


143.61 - 146.66 


258 


1592, 3184 


p-H 2 4 13 -3 22 , CO J=18-17, [Oi] 3 P - 3 Pi 


2 


B 


78.37 - 79.73 


147 


1592, 3184 


o-H 2 4 23 -3 12 , p-H 2 6 15 -5 24 , OH *IIi A i/a - 2 n 3/2 , 3/2 , CO J=33-32 


1 


R 


156.73 - 159.43 


239 


1592, 3184 


[Cn] 2 P 3/2 - 2 P 1/2 , p-H 2 3 3 i-4 04 


2 


B 


89.29 - 90.72 


121 


1592, 3184 


p-H 2 3 22 -2„, CH + J=4-3, CO J=29-28 


1 


R 


178.58 - 181.44 


204 


1592, 3184 


o-H 2 2 12 -l i, CH + J=2-l, o-H 2 2,,-2 12 



For the sources in our sample, this simple picture is compli- 
cated for two reasons: (i) in some observations the source is not 
centred in the central spaxel. In this case, even when observing 
an unresolved source, the emission is detected in a number of 
spaxels around the source position; (ii) the spectroscopic PSF 
width in the observed spectral ranges is equal or larger than the 
spaxel size (the PSF width for a source centred on the central 
spaxel vary between ~9" at 60 fim and ~13" at 180 fim), thus 
may cause line and continuum detection also in the outer spax- 
els. This effect is greater at the longer wavelengths, for bright 
sources, and when the source was not centred. 

A simple method to distinguish between extended and non- 
extended line emission is to calculate the line-to-continuum ratio 
in each spaxel of the grid. This is defined as the ratio between 
the line flux, integrated over its spectral profile, and the contin- 
uum, integrated over one spectral resolution element. If the line 
and the continuum emission originate from the same region, i.e. 
from the star-disk system, then the line and the continuum PSF 
will peak at the same position, and the line-to-continuum ratio 
is expected to be constant in all spaxels. If the line originates 
from a region more extended and/or offset with respect to the 
continuum emitting region, we should measure higher line-to- 
continuum ratios in the outer spaxels along the direction of the 
extended emission. Following this method, we compute the line- 
to-continuum ratios in all the spaxels for all the lines covered by 
our observations. In the [Oi] 63 /mi maps shown in Fig. [2] the 
presence of extended emission is indicated by line-to-continuum 
ratio which are larger than that measured on-source (i.e. where 
the continuum is at a maximum). 

To quantify and localise the extended emission identified 
through the line-to-continuum ratios, we subtracted on-source 
line and continuum emission, which could be detected in the 
outer spaxels because of the spectroscopic PSF width and/or off- 
centre sources. This analysis was applied to all the obtained line 
emission maps and allowed us to determine which lines show 
evidence for extended emission, to identify the spaxels where 
this is detected and at which confidence level. The mathematical 
formulation of the applied analysis is explained in Appendix iBl 
while the obtained continuum and residual line emission contour 
plots for the [0 1] 63 fim line are shown in Fig.[2](right panels). 



3. Results 

All of the observed sources show emission in a number of atomic 
([Oi] 63.2, 145.5 fim, [Cn] 157.7 fim) and molecular (high-J 
CO, ortho and para H 2 0, and OH) lines (see Fig.Q]and Tab.©. 
The spatial distribution of the atomic and molecular lines is dis- 
cussed in Sect. 13. H and Sect. 13.21 



3.1. Atomic [O \], [C ii] emission: correlation with optical jets 
and millimetre outflows 

Based on the analysis presented in Sect. I2.3I and Appendix [B] 
all of the sources in the sample show extended emission in the 
[Oi] 63 fim and in the [Cn] 158 fim lines, with the exception 
of RW Aur which shows faint and unresolved [C 11] emission. 
Despite the limitations imposed by the low spatial resolution 
of the PACS data, the extended emission detected in the bright 
[Oi] 63 fim line is spatially correlated with the direction of 
the jets/outflows. In the following we examine the [Oi] 63 fim 
maps obtained for each source in the sample, and compare 
these with observations at optical, near-infrared and millimetre 
wavelengths (see Fig.|2|i. 

TTau: 

T Tau is a multiple system consisting of the optically visi- 
ble northern component T Tau N, a nd the "infrared" c ompanion 
T Tau S located 077 to the south dDvcket al.lll982l) . The lat- 
ter is itself a close binary wi th a separation bet ween the com- 
ponents Sa and Sb of 0713 dKohler et al.l l2008). T Tau N has 
been classified as a Class II source, while T Tau S is deeply 
embedded and prob ably a Class I source dFurlan et al.1 [2006; 
iLuhman et all 1201 Oh . Both T Tau South and T Tau North drive 
jets whose forbidden optical line emission is detected up to a 
distance of ~40" along the North-South and the West-East di- 
rections (ISolf&Bohmt 19991) . 

PACS photometric observations at 70, 110 and 160 fim ac- 
quired within the GASPS project, and presented in Howard et 
al. (in preparation) show that the FIR continuum emission as- 
sociated with T Tau is extended (>4"). In the spectroscopic ob- 
servations presented in Fig.[2]the unresolved multiple system is 
centred on the central spaxel. 

The [0 1] 63 fim line and continuum emission are maximum 
on the central spaxel, and are detected in all the spaxels across 
the PACS field of view. The line-to-continuum ratios are larger 
in the outer spaxels, indicating that the line emission is more 
extended than the continuum. After subtracting the on-source 
line emission scaled to the detected continuum level in all 
spaxels (see Appendix [B] for details), we detected residual line 
emission above the 5<x confidence level in most of the outer 
spaxels up to a distance of ~28". The residual emission shows 
two peaks located -10" to the north-east (F res id ua i ~ 1.7 10 
Wnr 2 ) and ~8" to the south (¥ residual ~ 8.2 1(T 16 WtrT 2 ) with 
respect to the source. At the PACS resolution it was impossible 
to distinguish the emission associated with the two jets detected 
at optical wavelengths but only possible to identify the presence 
of extended line emission. Note that the bulk of the [0 1] 63 
fim line emission originates from a region of <11" around the 
source but extends up to ~30"(i.e. up to the edge of the PACS 
grid) similarly to what is seen in the optical forbidden lines 
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Table 3. Continuum emission, in Jy, in the spectral ranges covered by the Herschel/PACS observations. The continuum flux was 
estimated after integrating the emission over the 25 PACS spaxels. 



Source 


TTau 


DG Tau A 


DG Tau B 


FS Tau A+B 


RW Aur 


Spectral range (pm) 


Cont ±ACont (Jy) 


62.93 - 63.43 


132.4 


± 0.8 


14.6 ±0.3 


11.9 ±0.2 


5.2 ± 0.2 


0.7 ± 0.2 


71.82-73.33 


129.9 


±0.1 


15.5 ±0.1 


11.48 ±0.06 


4.65 ± 0.05 


2.25 ±0.10 


78.37 - 79.73 


130.5 


±0.2 


16.3 ± 0.2 


11.9 ±0.2 


4.4 ±0.1 


2.3 ±0.3 


89.29 - 90.72 


129.6 


± 0.2 


17.3 ± 0.2 


12.3 ±0.1 


5.14 ±0.10 


3.2 ±0.1 


143.61 - 146.66 


105.36 


± 0.10 


17.19 ±0.09 


16.93 ± 0.06 


6.69 ± 0.05 


1.36 ± 0.07 


156.73 - 159.43 


102.6 


±0.1 


17.3 ± 0.1 


17.43 ± 0.06 


7.19 ±0.05 


1.0 ±0.1 


178.58 - 181.44 


81.9 


t0.2 


15.4 ± 0.1 


15.57 ±0.09 


6.71 ±0.08 


0.8 ±0.1 



Table 4. Atomic and molecular line fluxes in W m 2 . The line fluxes were obtained by summing the emission over the spaxel where 
the line was detected (see text for details). 



Source 




TTau 


DG Tau A 




DG Tau B 




FS Tau A+B 


RW Aur 




Transition 










± AF(Wm- 


! ) 










[Oi] J P r j P 2 


63.184 


1.91 ±0.01 10" 14 


1.79 ±0.07 10- 15 


7.3 ± 0.4 IO" 


16 


5.2 ± 0.3 10" 


16 


2.0 ± 0.2 10" 


16 


[Ol] 3 P - 3 P! 


145.525 


8.8 ±0.3 10- 16 


8.6 ± 0.9 10- 


17 


2.4 ± 0.2 10- 


17 


2.3 ± 0.4 10" 


17 


1.3 ± 0.4 10" 


17 


[Cn] 2 P 3/2 - 2 P 1/2 


157.741 


7.5 ±0.3 10- 16 


3.0 ± 0.2 10- 


16 


2.4 ± 0.5 10" 


17 


4.2 ± 0.9 10" 


17 


<4 10- 18 




CO J = 36-35 


72.843 


1.0 ±0.1 io- 16 


1.5 ±0.4 10- 


17 


<5 10- |!i 




<4 IO-' 8 




< 1 io-" 




CO J = 33-32 


79.360 


2.1 ±0.2 IO" 16 


< i io-' 7 




<9 10-' 8 




<4 IO-' 8 




< i io- 17 




CO J = 29-28 


90.163 


3.9 ±0.2 IO" 16 


2.9 ± 0.9 10- 


17 


< 3 10-' 8 




< 3 IO" 18 




<5 IO" 18 




CO J = 18-17 


144.784 


1.29 ±0.03 IO" 15 


5.2 ± 0.9 10- 


17 


1.3 ±0.2 IO" 


17 


3.6 ± 0.4 10" 


17 


1.2 ± 0.2 10" 


17 


p-H 2 6 15 - 5 24 


78.928 


1.0 ±0.2 10-'" 


< 1 io-" 




<7 IO-' 8 




< 3 IO" 18 




< 1 10"" 




p-H 2 3 22 - 2n 


89.988 


3.3 ±0.2 10- 16 


< 1 io- 17 




< 3 IO" 18 




4.9 ± 1.6 10" 


18 


< 6 IO" 18 




p-H 2 4 13 - 3 22 


144.518 


7.0 ±2.9 10-' 7 


< i io- 17 




< 6 IO" 18 




< 3 10- 18 




<3 10- 18 




p-H 2 3 3 i - 4 04 


158.309 


< 1 io-' 7 


< i io- 17 




<6 10- 18 




<5 10- 18 




<4 10- 18 




o-H 2 8i 8 - 7 07 


63.324 


3.2 ± 1.2 10- 16 


<4 IO" 1 ' 




< 1 io-" 




<4 IO" 1 ' 




< 1 io- 1 ' 




o-H 2 7 07 - 6i6 


71.947 


2.9 ±0.2 IO" 16 


< i io- 17 




<5 10- 18 




<4 IO" 18 




<9 10- 18 




o-H 2 4 23 - 3i2 


78.741 


5.6 ±0.2 10^ 16 


1.9 ± 1.4 10- 


17 


<7 IO" 18 




1.2 ± 0.2 10" 


17 


2.1 ± 0.8 10" 


17 


o-H 2 2 12 - loi 


179.527 


6.1 ±0.2 IO" 16 


1.5 ±0.3 10- 


17 


2.6 ± 0.6 10- 


IS 


3.6 ± 0.4 10" 


17 


1.9 ± 0.6 10" 


17 


o-H 2 2 2 i - 2k 


180.488 


1.7 ±0.2 10^ 16 


< 1 io- 17 




1.9 ± 0.6 10- 


IS 


1.1 ± 0.5 10" 


17 


1.0 ± 0.7 10" 


17 


oh 2 n 1/2il/2 


79.110 


8.9 ±0.1 lO" 16 


5.2 ± 0.5 10- 


17 


8.7 ± 2.4 10- 


18 


1.8 ± 0.1 10" 


17 


7.7 ± 3.2 10" 


18 


oh 2 n 3 / 2 3 / 2 


79.180 


1.18 ±0.01 IO" 15 


2.8 ± 0.5 10- 


17 


6.0 ± 2.4 10" 


IS 


1.7 ± 0.1 10" 


17 


1.7 ± 0.3 10" 


17 



(ISolf & Bohmir i999). 
DG Tau A: 

DG Tau A is a strongly accreting Class II source (e.g., 
lHartigan et al.|[l995l iLuh man et al"] |2010l) associated with a jet. 
The DG Tau A jet wa s first detected at optical wavelengths by 
Mundt&FrielJllIl. High angular resolution studies at optical 
and NIR wavelengths showed a bright collimated blue-shifted 
lobe moving with radial velocities up to ~350 kms" 1 , whilst 
only faint emission was detected on t he red-shifted side (e.g., 
iDougados et"ai1l2000t IPvo et alj 12003). Spectroscopic observa- 
tions taken with the H ub ble Space Telescope and presented by 
IBacciotti etal] (120001) and lMaurri et all (l2012h indicated that the 
bulk of the emission in the [O i] 6300 A, and in the other optical 
forbidden lines ([S n] and [Nn] lines) comes from the first 1'.'3 
of the blueshifted jet lobe. Beyond 1'.'3, two strong bow-shocks 
are detected. One is at ~3"-4'.'5 from the source and the other is 
at ~9-10" from the source. Both working surfaces are moving 
with a proper motion of ~0.3"/year. 

Fig. |2] (left panel) shows that in our PACS observations the 
source is not centred in the central spaxel but lies 6'.'7 to the 
east with respect to it. Thus, both the line and the continuum 
emission are detected in a number of spaxels around the source 
position. However, whilsty the continuum emission peaks at 
the position of spaxel (1,1), the [Oi] 63 fim line peaks at the 



position of spaxel (2,2), i.e. ~6'.'7 to the west with respect to 
the source. This indicates that the line emission is offset with 
respect to the continuum, and is possibly more extended. The 
right panel of Fig. [2] shows the residual line emission offset 
with respect to the continuum emission, and displaced along the 
direction of the blue lobe that is detected at optical wavelengths. 
The residual line emission reaches its maximum ~1" from the 
source (¥ res iduai ~1.5 IO -15 WirT 2 ) and extends ~11". 

DG Tau B: 

DG Tau B has b een classifies as a Class I source by 
iLuhmanet ail d201Q|) ; iRebull et all d2010t) and is associated 
with the bright, asymme tric HH 159 iet dMundt & Friedlll983b 
lEisloffel & Mundtl(l998l) . The red lobe consists of a chain of 
bright knots detected in the optical forbidden lines ([Nn], [Oi], 
[S n]) that extends ~55" from the source. The fainter blue lobe 
is only detected up to ~10" from the source in the same lines. 
Observations at millimetre wavelengths in the CO low-J lines 
confirm the asymmetric structure showing a slow, wide-angle 
red-shifted outflow which is displaced along the axis of the op- 
tical jet, w hile blue-shifted CO emis sion is faint and confined 
on-source (Mitch ell et alll 1994 1 19971) . The non-detection of the 
blue lobe at millimetre wavelengths suggests that the observed 
asymmetry is not due to obscuration of the blue lobe but to dif- 
ferent physical conditions in the two lobes, where the bright red 



5 



L. Podio et al.: HerschelfPACS observations of young sources in Taurus: the far-infrared counterpart of optical jets 




Wavelength (\im) 

Fig. 1. Spectral ranges covered by the HerschelfPACS observations for all of the sources in the sample. The wavelengths of the 
targeted atomic ([O i], [C n]) and molecular (H2O, OH, CO) lines are shown by the vertical dotted red lines. 



lobe is two times slower and less i onized but denser , and more 
collimated than the faint blue lobe (Podio et al. 2 01 lh . 



The [Oi] 63 fim emission map in the left panel of Fig. [2] 
further confirms this asymmetric structure. The peak of the line 
emission (in spaxel (2,2)) is not coincident with the peak of the 
continuum emission (at spaxel (1,1)), and the line emission is 
displaced ~34" (i.e. up to the edge of the PACS grid) along 
the direction of the optical-jet/mm-outflow red lobe. Residual 
emission above the 5<x level was detected in four spaxels along 
the red-lobe jet PA, decreasing from 1.5 10~ 16 WrrT 2 at ~1" 
down to 2 10~ 17 W irT 2 at 26" from the source (see Fig. [2] right 
panel). An order of magnitude decrease in flux between the first 
few arcseconds of the jet and the emission at 20 "-30" from 
sourc e was also seen in the optical forbidden lines ( Podi o~et~a"fl 
1201 lb . Therefore, the PACS maps suggest that the [Oi] 63 /im 
line is tracing a warm outflow component, intermediate between 
the hot, fast and collimated jet traced by the optical forbidden 
lines and the cold and slow wide-angle outflow traced by 
millimetre CO emission. Higher angular and spectral resolution 
observations are required to fully analyse the spatio-kinematical 
structure of this warm component. 



FS Tau A, B: 

FS Tau A is a Class II clo se binary system (separation 
~(y/23 bv IWhite & Ghezl I200H class by iFurlan et al.1 120061 
iLuhman et al.l 120101) . associated with a bright reflection neb- 
ulosity but showing no clear evidence of an associated jet 
(lKrist etal.lll99i lEisloffel & Mundl T998). FS Tau B is an em- 
bedded not o ptical ly visible Class I source dLuhman et al.ll20l61; 
iRebull et al.l 12010)" located ~20" to the west. This source is 
driving a parsec scale c ollimated outflow (M undt et al.l fl984; 
lEisloffel & Mundtlll998l) . 

The [Oi] 63 fim map in the left panel of Fig. [2] shows 
two continuum peaks, at the position of spaxel (2,2) and (1,4) 
corresponding to the position of FS Tau A and FS Tau B, 
respectively. Line emission was detected in a number of spaxels 
located al ong the PA of the o ptical jet associated with FS Tau B 
(PA~55°, Mu ndt et all 1 1984 and around FS Tau A. However, 
because of the limited spatial resolution of PACS, it is not 
possible to separate the line and continuum emission associated 
with the two sources and to apply the analysis presented in 
Appendix iBl Thus, the total line and continuum emission from 
the FS Tau system (Aa+Ab+B) is considered throughout the 
paper. 
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T Tau: [SII] 6731 A,+ [01] 63.2iim h T Tau: [01] 63.2UIT1,, continuum + line residual 
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Fig. 2. Left panel: [O i] 63 fim line spectra in the 25 spaxels of the PACS array. Each spaxel is labelled by its (x,y) index and its 
position is indicated by a black cross. The RA and Declination offsets with respect to the source (in arcseconds) are indicated on the 
x and y axis. The black star is the source position and the blue/red dashed lines the position angle of the optical-jet blue/red lobe. In 
each spectrum, vertical and horizontal dotted lines show the line wavelength and the zero flux level, respectively. In those spaxels 
where the line is detected above the 3<x level, the spectrum is drawn with a black thick line, and the estimated continuum level (in Jy), 
line flux (in W rrT 2 ), and line-to-continuum ratio (in red) are indicated. Where the line is not detected only the estimated continuum 
level (in Jy) is given. The opti cal-jet e mission in the rSn~|/bl671 6, 6731 A forbidden lines is overplotted (gr een contours). Contours 
are from lSolf & Bohml (fl999h (T Taul. lEisloffel & Mundtl (fl998h (DG Tau A, DG Tau B, FS Tau A+B), and lDougadosltal] {2000) 
(RW Aur). Right panel: Contour levels of continuum (dotted lines) and line residual (solid lines) emission obtained from the analysis 
in Appendix IB"1 The spaxels where residual line emission is detected with a confidence level >5 are highlighted by a red circle (for 
RW Aur the red circles indicate residual emission above the 2cr level). The contours indicate that the line emission is shifted and/or 
more extended with respect to the continuum emission along the optical jet PA. 



RW Aur: 

RW Aur is a binary syste m (~1'.'4 separation) and has been 
classi fied as a Class II source (I White & Ghezll200UlFurlan et al.l 
I2006h . RW Aur A is associated with an asymmetric jet whose 
bright red lobe is detected in the optical [Sn], [Nn l, and TOil 
forbidden lines up to ~ 15" from the source (e.g.. |Hirth et all 
1 19941 119971: iDougados et all2000t iMelnikov et alJ2009h . 



The binary system is not resolved with PACS. The emission 
map at 63 fim in Fig. [2] shows that both the continuum and line 
emission peak on the central spaxel but the [O i] 63 /vm line is 
detected also in a few outer spaxels lo cated along the di r ection 
of the optical jet red lobe (PA~120°, iHirth et al.l fl994l Il997t 
IDougados et al.l 120001) . This suggests the line is tracing the 
warm FIR counterpart of the hot optical jet. After subtracting 
on-source emission there is no residual emission above 5cr. 
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However, residual emission above 2<x is observed along the 
red-shifted jet direction at ~9" and ~15" from the source 

{^residual ~L8 lO^'Win- 2 ). 

Since all of the sources in our sample show extended 
[Oi] 63 pm emission spatially correlated with the direction 
of the optical jets we also expect the line profile to be blue- 
or red- shifted in agreement with what is observed at optical 
wavelengths. However, the low spectral resolution offered by 
PACS (-88 kms 1 at 63 pm) is further hmited by the fact 
that extended emission can broaden, and shift, the line spectral 
profile (see PACS manual). In the case of T Tau, DG Tau A, and 
FS Tau A+B this effect clearly dominates, as we can see that 
the line peak shifts from red to blue velocities moving across 
the PACS field of view (see Fig. [2j. However, in the case of DG 
Tau B and RW Aur the [Oi] 63 pm spectral profile is much 
larger than the instrumental one (FWHM up to 200 kms 1 ) and 
consistently shows a red-shifted peak velocity in all the spaxels 
displaced along the direction of the optical jet red lobe. The 
peak velocity estimated from those spaxels showing residual 
line emission >5<x is of ~45 kms 1 and ~95 kms 1 for DG 
Tau B and RW Aur, i n agreement with values estimated from 
optical forbidden lines dPodio et al.ll201 lblMelnikov et al.ll2009h . 

Even if for all of the sources in the sample extended emis- 
sion is detected in the atomic [O i] and [C n] lines suggesting a 
jet/outflow origin, the emission from the disk may dominate on- 
source. For this reason we compare observed atomic lines with 
predictions from both disk and shock models in Sect. [4] 

3.2. Molecular H 2 0, CO, and OH lines: unresolved emission 
and high excitation lines 

In our brightest target, T Tau, we have detected five ortho-H20 
and four para-H20 lines, including lines from highly excited 
levels, such as the 0-H2O818 - loi line (E„ p ~1070 K). In the 
other sources, only a few of the lower excitation H2O lines were 
detected, up to the 0-H2O423 - 3 12 line (E up ~432 K). T Tau and 
DG Tau also show emission from high excitation CO levels up to 
CO J=36-35 line (E llp -3668 K). Only the lowest CO transition 
covered by our observations, i.e. the CO J=18 
K), is detected from our other sources (see Tabled 

We have also detected the J=5-4 CH + emission line at 72.14 
pm from T Tau. The other CH + transitions falling in the spectral 
range covered by our observations (i.e. the J=4-3 and J=2-l lines 
at 90.01 and 179.60 pm, respectively) are blended with the p- 
H 2 and o-H z O lines at 89.988 and 179.527 pm. These are the 
brightest transitions of para and ortho water in all of our sources. 
Because of the low spectral resolution of PACS and that the wa- 
ter lines are an order of magnitude brighter than the detected 
CH + emission line (~4 10~ 17 W irT 2 ), we cannot deconvolve the 
lines and recover the flux of the faint CH + lines. The origin of 
the spectrally and spatially unreso lved CH + line i s unclear. This 
line could be excited in the disk (iThi et alJl201 lh and/or in the 
directly irradiated outflow cavity walls (Bruderer et alJl2010l) . 

On the contrary to that found for the atomic [0 1] and [C 11] 
lines, the emission in the molecular lines is both spectrally and 
spatially unresolved with PACS. By applying the analysis in 
Appendix iBl we found that, even when line emission is detected 
out of the central spaxel, the line-to-continuum ratio is constant 
across the PACS grid and there is no residual line emission above 
5cr after the subtraction of on-source line and continuum. Hence, 
the line and continuum emission detected in the outer spaxels are 
due to the spectroscopic PSF and the line and continuum emit- 
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Fig. 3. The fluxes of the [Oi] 63 /mi, the [Cn] 158 /mi, the 
CO J=18-17, and the 0-H2O 179.5 /mi lines are plotted versus 
the source luminosity {left panel) and their adjacent continuum 
(right panel) for all the sources in our sample (black, red, green, 
and blue dots). For FS Tau A+B and T Tau N+S the plotted 
source luminosity is a lower limit, since L* of the Class I com- 
ponents FS Tau B and T Tau S is unknown. The upper limit on 
the stellar luminosity for D G Tau B corresponds to its bolometric 
luminosity as estimated by Kruger Tt alJ (1201 lh . 



ting regions are unresolved. This implies that the atomic and 
molecular lines have different spatial distributions and origins. 
While the extended atomic emission is clearly associated with 
the optical jets, the compact molecular lines may be excited in 
the disk, in the UV-heated outflow cavities, and/or in the first 9" 
(~1300 AU) along the jet. To understand which of these compo- 
nents dominates the line emission, we investigate possible sce- 
narios in the following section. 
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Fig. 4. Observed [O i], [C n], CO, and H2O line fluxes and ratios (re d dots) are compared with predictions from a subsample of disk 
models from the DENT grid (Wo itke et al.ll2010llKamp et al.ll201 ll) . Contours encircle 85% of the DENT disks for star luminosity 
values from 10~ 2 L to 9 L Q (red, yellow, green, cyan, blue, and black crosses). The total errors on the observed values (dotted 
red lines) are obtained summing the error due to the line signal-to-noise (solid red lines) and the 30% error affecting PACS flux 
calibration. The subsample of disk models is obtained by selecting M<2.5 M , T e ff < 5500 K, L<9 L Q , R,„=R., 1( ;,/, and dust-to-gas 
ratio = 0.01. 



4. Discussion 

Fig. [3] shows the fluxes of the brightest observed lines 
([Oi] 63.2 ^m, [Cn] 157.7 pm, H 2 179.5 pm, CO J=18-17) 
versus the source luminosity (left panel) and the continuum flux 
adjacent to the considered line (right panel). The plots show no 
clear correlation, except for the [Oi] 63 pm line. Howard et 
al. (in preparation) examine a large sample of Taurus sources, 
including optical-jet sources (i.e. sources with jet signatures, 
such as forbidden emission lines which are extended and/or 
show blue-shifted profiles), and sources showing no evidence 
of outflowing activity. Most of the Class II jet-sources in Taurus 
are associated with micro-jets extending only up to a few arc- 
seconds away from source in the typical optical tracers. These 
sources show unresolved [Oi] 63 pm emission at PACS reso- 
lution, making it difficult to disentangle jet and disk emission. 
However, Howard et al. (in preparation) find that optical-jet 
sources show excess [Oi] 63 pm emission with respect to the 
tight F ([Oi] 63 ^m) - Continuum (63 pm) correlation which is 
found for sources with no evidence of outflowing activity. By 
using this correlation and the measured continuum fluxes we 



find that the emission from the disk may account for 3% to 15% 
of the observed [Oi] 63 pm line flux. This suggests that in jet 
sources emission from the jet/outflow dominates over possible 
disk emission in the [Oi] 63 pm line. 



To verify this hypothesis and constrain the origin of the unre- 
solved molecular emission in the following sections we compare 
the atomic and molecular line fluxes, and their ratios, with disk 
models and shock models predictions (Sect. l4,ll and Sect. 14. 21 re- 
spectively). The line ratios are computed by using the line fluxes 
summarised in Table [4] Unfortunately, an analysis of the varia- 
tion of line ratios with distance from source is not possible from 
several reasons, i.e. (i) in some observations the source is not 
centred on the central spaxel; (ii) because of the spectroscopic 
PSF the emission in the outer spaxels is strongly contaminated 
by on-source emission overlapping to local extended emission; 
(iii) the source position on the PACS field of view may be differ- 
ent when observing at different wavelengths. 
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4.1. Emission from disks 

To understand how much the disk can contribute to the extended 
[Oi] and [Cn] atomic lines and to the unresolved molecular 
emission, we compare observed line fluxes and ratios with pre- 
dictions from the DENT grid of dis k models dWoitke et al.l2 010; 
iPinte et alJlIoToHKamp et al.ll201 ll) . The DENT grid consists of 
300 000 disk models spanning a large range of parameters defin- 
ing the source (mass, M», temperature, T e ff, luminosity, L», UV 
excess) and the disk (disk gas mass, gas-to-dust ratio, inner and 
outer disk radius, R,„, R OM ,, surface density, flaring, dust grain 
size distribution, dust settling, disk inclination) properties. For 
each model, the dust temperatures, gas temperatures and chemi- 
cal structure are computed to produce a large set of observables 
such as the SED and selected FIR and submillimetre line fluxes. 
We consider a sub-sample of models in the grid corresponding 
to typical low-mass young stellar objects (YSO) and T Tauri star 
properties (M, < 2.5 M , T eff < 5500 K, L» < 9 L ). We also 
assume that the inner disk radius, R,„ , is located at th e position 
of the dust sublimation radius, R su bi dPinte et al.ll2008l) and that 
the dust-to-gas ratio is 0.01, as expected in the case of young 
primordial disks. 

Fig. 3] indicates that the selected disk models cover a large 
range of line ratios but cannot reproduce the absolute fluxes of 
the observed bright lines. In particular, we observe [Oi] 63 fim 
fluxes up to 1.9 10~ 14 WirT 2 and 1.8 10~ 15 WirT 2 for T Tau 
and DG Tau A, which cannot be produced by the disk even for 
high star luminosity. DG Tau B and FS Tau A+B show lower 
[Oi] fluxes (-5-7 10~ 16 WirT 2 ) but, if we consider that their 
luminosity is lower than 1 L , disk models cannot account for 
the observed flux. The emission in the fine structure lines from 
disk models is similar to that predicted by PDR models, with the 
fundamental difference that it originates over a wider range of 
density/temperature. The [O i] lines are optically thick and close 
to LTE (the critical densities for the [O i] 145 um and [O i] 6 3 fim 
lines are 6 10 4 , 5 10 5 cm -3 respectively: iKamp et alj|2010l) . 

The discrepancy between observed fluxes and predictions 
from disk models is even more evident for the molecular lines: 
while observed CO 18-17 and H2O 179.5 fim line fluxes are of 
10 17 -10 15 WirT 2 and 3 10~ 18 -6 10~ 16 WirT 2 , the predicted 
fluxes are always lower than a few 10~ 17 WirT 2 for the CO 18- 
17 line and 10~ l7 Wm" 2 for the H 2 179.5 fim line. According 
to the disk models these lines originate at the disk surface and 
the line excitation temperature determines the radial extent of 
the emitting area. For solar-t ype stars these wa ter lines are op- 
tically thick and not in LTE dAresu et al.ll2012l) . The high-J CO 
lines are generally in LTE due to their low critical densities, and 
they can be optically thick depending on the details of the model. 
Higher atomic and molecular line fluxes can be obtained by as- 
suming a lower dust-to-gas ratio (dust-to-gas = 0.001) and very 
large inner radii (e.g., R,„ >10 R sm m). These values, however, 
are typical of more evolved disks but are not appropriate to de- 
scribe the young disks observed around Class I and II sources. 
Since some of the observed sources are strong X-rays emitters 
dGiidel et al.ll2007l) . we also checked the effect of X-rays o n the 
considered emission lines in the work of Aresu et al .1 (1201 lh and 
Meii erink et alj d2012l) ; lAresu et all d2012ir It can be shown that 
for the considered subsample of disk models the line fluxes are 
not significantly enhanced by X-rays (less than a factor of 2). 



4.2. Emission from shocks 

The strong spatial correlation between the [0 1] 63 fim line and 
the optical jets and the large observed fluxes suggest that the 



o 0.1000- 




[OI] 63.2 \im/ [OI] 145.5 

Fig. 5. Observed atomic line ratios are compared w ith predic- 
tions from PDR (dotte d lines. iKaufman et aljll999h and fast J- 
type shock (solid lines. iHollenbach & M cKee 1989) models. For 
the PDR models, the labels on the broken curves indicate the gas 
density (n) and the intensity of the FUV field (Go), respectively. 
For the shock models, the pre-shock density is marked at the 
lowest shock-velocity point (V s h oc k=30 kms _I ), and the shock 
velocity increases along the full curves, up to V s hock= 150 km s _1 . 



emission in the atomic lines originates from the shocks occurring 
along the jet and/or from the UV-heated gas in the outflow cavity 
walls, rather than in the circumstellar disk. Thus, we compare 
atomic line fluxes and ratios with predictions from shock models 
and PDR models. 

Fig. [5] shows the observed atomic line ra- 
tios [Oi]63yum/145yum, hereafter [Oi]63/145, and 
[Cn]158^m/[Oi]63^m, hereafter \CiiV\6i\ and the pre- 
dictions by PDR models (Kaufm an et al.lfl9 99) and fast J-type 
shock models, in which H2 is fully diss ociated and the gas is 
partially ionized in a radiative precursor dHollenbach & M cKee 
11989b . The lowe r-velocity (V^-lO^Okm s " 1 ) C- and J-type 
shock models of iFlower & Pineau Pes Foretsl (|2010), which do 
not incorporate a radiative precursor, are absent from this plot, 
as they predict [Cn] 158 fim emission from 2 to 8 orders of 
magnitude lower than the [0 1] 63 fim emission with increasing 
pre-shock density. 

Fig. [5] indicates that both PDR models with densities > 
10 4 crrT 3 and FUV field Go > 10 3 and dissociative J-type shocks 
with low pre-shock density (n~10 3 -10 4 crrT 3 , V shock = 30-150 
km s _1 ) are able to reproduce the [O i]63/145 and the [C n]/[0 1] 
line ratios for most of the sources in our sample. The exceptions 
are RW Aur, which cannot be reproduced by PDR models be- 
cause [Cii]/[Oi] <0.02, and DG Tau A, which cannot be repro- 
duced by the fast J-shock models because of its large [C n]/[0 1] 
line ratio (-0.17). 

The [Oi] 63 fim maps presented in Fig. [2] show that the 
[0 1] line emission originates from a region extending up to 30" 
from the source. Such extended emission could not be produced 
by UV-illuminated outflow cavities. Thus a shock origin likely 
dominates the [0 1] line emission. However, part of the observed 
emission may originate from the illuminated outflow cavities, 
which would also explain [C n]/[0 1] line ratios higher than what 
is predicted by shock models. Note that the only object showing 
no evidence of surrounding cloud material or outflow cavities is 
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o-H 2 180.5 fim / o-H z 179.5 fim CO J = 29-28 / CO J = 1 8- 1 7 

Fig. 6. Observed H2 O and CO line ratios (black dots) ar e compared with the predictions of slow C-type (red lines) and J-type {green 
lines) shock models (Flo wer & Pineau Pes Foretsll2010h . The solid lines indicate the errors due to the line signal-to-noise, while the 
dashed lines indicate the total error, obtained by adding the 30% error in the PACS flux calibration. The crosses along the red and 
green curves correspond to increasing shock velocity, from 10 km s _1 to 40 km s _1 (C-type shocks) or from 10 km s _1 to 30 km s _1 
(J-type shocks). At the lowest velocity (V s h oc k=10 km s _1 ), the pre-shock density is given. 



p-H^O lines o-H O lines P"H,0 lines o-H ; lines 




78.9 90.0 144.5 158.3 63.3 71.9 78.7 179.5 180.5 78.9 90.0 144.5 158.3 63.3 71.9 78.7 179.5 180.5 

Line Wavelength ([im) - Arbitrary scale Line Wavelength (fim) - Arbitrary scale 



Fig. 7. Observed H2O line r atios for T Tau (red points/arrows) are compared with the predictions of C-type {left panel) and J-type 
(right panel) shock models dFlower & Pineau Pes For ets 2010) with relatively high pre-shock densities . The H2O line ratios and 
upper limits for the other sources in the sample are overplotted (black points and arrows). 



RW Aur, which is also the only source in our sample with no 
detected [Cn] 158 pm emission. 

Petermining the origin of the observed molecular lines is 
more difficult because, as explained in Sect. 13.21 they are spec- 
trally and spatially unresolved with PACS. The fact that the large 
CO and H2O line fluxes cannot be reproduced by disk mod- 
els, with parameters typical of low mass YSO and T Tauri stars, 
favours either a shock origin or a PPR origin in UV-heated out- 
flow cavities. The molecular emission is more compact than the 
atomic emission but could, nonetheless, arise in a shock or an 
outflow. 

Evidence of molecular emission associated with a jet, but 
which is less extended than the atomic emission, has been found 
in evolve d CTTSs, o b served at near-infrared wavelengths. In 
particular, iBeck et a"D (|2008) analysed high angular resolution 
observations (~(J.'l) of the H2 2.12 pm line in CTTSs that are 
driving jets. They showed that the emission in this line is more 
compact than the emission in atomic optical forbidden lines 
([Oi], [S n], [Nil]) but still exten ds up to l" - 2" an d is spatially 
associated with the jet direction. IBeck et ail d2008l) investigated 
the origin of the H2 emission in three of the sources analysed 



in this paper (T Tau, PG Tau A, and RW Aur). They detected 
molecular hydrogen at distances >50 AU from the star and de- 
rived excitation temperatures >1800 K. They also found that 
the emission lines were consistent with existing shock models. 
Based on these measurements and the kinematics of the features, 
they concluded that most of the H2 toward these stars arose from 
shocks associated with the known HH objects rather than from 
quiescent disk gas illuminated by the central star. Moreover, a 
few recent studies have shown that high excitation H2O lines 
and even higher-/ CO transitions (up to J = 46-45) can be 
produce d in the outflows emana t ing from young Cla ss and I 
sources (Ivan Kem pen et ail l2010l;lHerczeg et al. 2012). 

In view of the results above, we have attempted to simulate 
the observed molecular line emission by me ans of shock models. 
As shown by iHollenbach & McKeel (1 19891) . fast J-type shocks 
give rise to a radiative precursor and are strongly dissociative; 
for a given [0 1] 63 pm line flux, the emission in the H2O lines 
is negligible. Accordingly, we have compared the observed H2O 
and CO line ratios with the predictions of slow C-type (V shock = 
10—40 kms' 1 ) and J-type (V shock = 10-30 kms -1 ) shock mod- 
els of iFlower & Pineau Pes Foretsl (l2010h (see Fig. [6] and |7). 
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In these shock waves, much of the mechanical energy is trans- 
form ed into H?Q and CO line radiatio n, as illustrated in Table 
1 of iFlower & Pineau Pes Foretsl (2010). Figs. [6] and [7] suggest 
that relatively high pre-shock densities are required to reproduce 
the intensities of transitions from high excitation levels of H2O , 
such as the o-H 2 8i 8 - 7 07 line (E up -1070 K). 

In the case of T Tau, all the H2O lines falling in the spec- 
tral range covered by our observations are detected with a good 
signal-to-noise ratio. A satisfactory fit to the observed line ratios 
is obtained for a C-type shock with a pre-shock density of 10 6 
cm" 3 and a velocity in the range 20-30 km s" 1 (see Fig . |7J) . The 
diameter (~ 220-360 AU) of the emitting region t hat we estimate 
for T Tau is comparable with that deduced by Spinoglio et alJ 
(2000), using a large velocity gradient (LVG) model to fit the 
observed FIR lines. 

For the other sources, we detected a maximum of four H2O 
lines. Under these circumstances, it makes little sense to look 
for a best-fit model; but the few line ratios and the upper lim- 
its overplotted in Fig.|7]indicate that C- and J-type shocks with 
pre-shock densities of 10 5 -10 7 cm" 3 and 1 4 — 1 5 cirT 3 , respec- 
tively, produce H2O line ratios which are consistent with those 
observed. Moreover, the models can easily reproduce the large 
H2O 179.5 fim line fluxes for an emitting region with a diame- 
ter of a few tens to a few hundreds of AUs (see Tab.[5]l, which is 
consistent with the fact that the source of the molecular emission 
is compact and unresolved with PACS. 

The observed 0-H2O 179.5 fim/CO 18-17 line ratios are 
lower, by up to a factor 4, than those predicted by the slow 
shock models: H 2 0/CO obs ~ 0.2-2; H 2 0/CO sh ocks > 0.8 (see 
right-hand panel of Fig. [6). Furthermore, the observed OH 
line fluxes are much larger (by up to a factor 10) than those 
predicted by these same models. H2O abundances which are 
lower than predicted by slow shock models hav e been found by 
Santangelo et al. (2012) and IVasta et alJ (1201 2h . from the anal- 
ysis of a large number of H2O emission lines associated with 
outflows in Class sources. These authors propose that the wa- 
ter ab undance may be reduced by UV dissociation (Bergi net alJ 
1 19981) and/or depletion onto grains. Photodissociation of H2O to 
OH, by the stellar FUV radiation field, is a po ssible explana- 
tion of these observations, as suggested also by Spinogli o et alJ 
(2000); but the further photodissociation of OH to O should also 
be taken into account. For further details, see Sect. 14.31 

From the comparisons with the models, we conclude that 
the atomic and the molecular emission - which have differ- 
ent spatial distributions (extended versus compact) - arise in 
shock waves with different characteristics (C- or J-type, with 
or without a radiative precursor). This conclusion is consistent 
with the results of previous analyses of FIR ISO observations 
of Class and I sources, which demonstrated that a single gas 
component cannot reproduce both the atomic and the molec- 
ular emission. It has been suggested that a dissociative J-type 
shock, occurring at the apex of the jet, and non-dissociative 
C-type shocks, occurring in the wings of the bow-like flow, 
may be responsible for exciting the atomic and the molecular 
lines, respectively (Nisi ni et al I 12001 . Alternatively, the emis- 
sion might arise from the UV-heated gas in the outflow cavity 
walls and small-scale C-type shocks occurring along the cavi- 
ties (van Kempen et al. 2010). We cannot exclude the possibility 
that the disk contributes to the molecular line emission, particu- 
larly t o the higher excitation water line s, as suggested by a recent 
study dRiviere-Marichalar et"aT]|2012l) . 

Follow-up observations with HIFI and ALMA are planned, 
with the goal of resolving the molecular lines, spectrally (with 
HIFI) and spatially (with ALMA), and identifying the contribu- 



Table 5. Diameter of the H2O emitting area (in AU) estimated 
from C-type shock models. 



source 


R(ITO) 




(AU) 


T Tau 


220 - 360 


DG Tau A 


35 - 270 


DG Tau B 


14-110 


FS Tau A+B 


54-412 


RW Aur 


39 -298 



tions to the emission from the disk, UV-heated outflow cavity 
walls, and shocks. 

4.3. T Tau: comparison with previous ISO data 

T Tau is the only source in our sample previously observed with 
ISO. The ISO observations were acquired with both the Short 
Wavelength Spectrometer (SWS) and the Long Wavelength 
Spectrometer (LWS) providing complete spectral coverage from 
2 fim to 190 fim. These observations show ed several H2O, CO 
(up to J„„=25), and OH em ission lines dvan den Ancker et alJ 
119991; ISpinoglio et"al]|2000l) . Our HerschelfPACS observations 
complement the previous ISO dataset by revealing high-J CO 
lines (up to J=36-35) and H2O lines which were not detected 
with ISO because of its lower sensitivity. 

The absolute line fluxes measured by ISO-LWS and 
Herschel are in very good agreement as shown in Fig. [8] 
Spinoglio et alJ d2000h modeled the observed lines by using an 
LVG code in a plane parallel geometry and found that most of 
the observed molecular emission may be explained by a dense 
and warm gas component (n=10 5 -10 6 cm" 3 , T=300-900 K) with 
an emitting area of diameter of 300-400 AU. 

Following these results we have modelled the full 
ISO+Herschel dataset with the shock models presented in the 
previous section. The left panel of Fig. [8] shows that most of 
the water lines are well fit with a single C-shock model (n=10 6 
cm" 3 , V J / !OC < : =20-30 kms" 1 , diameter of the emitting area of 
~360 AU). The exceptions are for a few high excitation water 
lines between 25 and 50 /mi observed with the ISO-SWS. To fit 
these lines, the pre-shock density must be about an order of mag- 
nitude larger (n=10 7 crrr 3 ). However, this discrepancy between 
models and observations may be partially due to intercalibration 
problems between the SWS and LWS (the continuum flux mea- 
sured by SWS is -10% higher than in th e LWS). 

The CO lines detected with ISO by ISpinoglio et"al] d2000h 
indicated that a warm component and a hot component are re- 
quired to reproduce all of the CO lines (Fig. [8]). This is further 
confirmed by our Herschel/PACS observations of CO lines up to 
J=36-35. While the very high-J CO lines detected by Herschel 
(CO J=36-35, 33-32, 29-28) are well fitted by the same C-shock 
model used to reproduce the H2O lines (we assumed that the 
H2O/CO abundance ratio is lowered by ~0.55 due to FUV ir- 
radiation of the shocked region) a warm gas component is re- 
quired to reproduce the lower- J CO lines (down to J, (/ ,= 14). 
We tentatively fit all CO lines by adding a warm gas com- 
pon ent using a RADEX LV G model in plane parallel geome- 
try dvan der Tak et alj|2007|) and choosing param eters which are 
similar to those used by Spinoglio et al. (2000) (n=10 6 cm" 3 , 
T=300 K, A=360 AU, N=1.5 10 18 cm" 2 ). This warm gas com- 
ponent also reproduces the OH line fluxes for an OH column 
density of ~5 10 16 cm" 2 . 
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Fig. 8. All the available ISO (black points) and Herschel (red points) observations of H2O {left panel) and CO {right panel) lines 
for T Tau are shown. Left panel: Line fluxes predicted by a slow C-shock for an emitting area of ~360 AU diameter and pre-shock 
density of n=10 6 cm" 3 (velocity V=20-30kms" 1 , dotted lines), or n=10 7 cm -3 (V= 10-20 km s" 1 , dotted lines). The low-excitation 
water lines are well reproduced by C-shock with pre-shock density of 10 6 cm" 3 while the higher excitation lines require higher pre- 
shock density. Right panel: CO line fluxes predicted by: a C-shock model (n=10 6 cirT 3 , V=20-30 kms -1 , emitting area diameter, 
D=360 AU, red line); an LVG warm gas component (n=10 6 , T=300 K, D=360 AU, N co =1.5 10 18 cm" 2 , blue line); the sum of the 
two models {green line). While the high-J CO lines are well reproduced by the same C-shock model used for water lines, a warm 
gas component is required to simultaneously fit the lower-J lines observed with ISO (down to J„ p =14). 



More detailed modelling of all the detected lines could 
be made by considering UV-heated gas in the outflow cavity 
walls, and small scale C -shocks along them, as discussed by 
Ivan Kempen et al.1 d2010l) . However, such detailed modelling is 
beyond the scope of the present paper. 

4.4. FIR cooling: an evolutionary picture 

To evaluate the efficiency of the FIR outflow component asso- 
ciated with the jet sources in our sample we estimate the total 
luminosity radiated away in the far-infrared lines, L (FIR). 

The [Oi], [Cn], and OH luminosities were estimated from 
the observed line fluxes (L [0 1]= L ([0 1] 63 + 145 pm), L [C n] 
= L [Cn] 157 pm, L (OH) = L (OH 79.11 + 79.18 pm)) 
while for the CO and H2O luminosity we give a lower limit, 
inferred by summing the fluxes of the observed lines, and an 
upper limit inferred by adding up the predicted line intensi- 
ties for all the transitions considered in the shock models (i.e. 
45 levels of 0-H7O, 45 levels of r >-H 7 Q, and 41 levels of CO, 
iFlower & Pineau Pes Foretsl l2010h . In particular, we consider 
the highest luminosity value obtained from the different shock 
models which reproduce the observed line ratios (i.e., J-shocks 
with pre-shock densities of no=10 4 -10 5 cm" 3 and C-shocks with 
no=10 5 -10 7 cm" 3 ). Then the cooling in the molecular lines is 
computed as L mo , = L (OH) + L (H 2 0) + L (CO), while the 
total cooling in FIR lines is estimated by summing up the cool- 
ing in all the detected species: L (FIR) = L [Oi] + L (OH) + 
L (H 2 Q) + L (CO). Note that followi ng previous ISO studies 
dGiannini et al.l200ltlNisini et alj2002}) we neglect the luminos- 
ity of the [C 11] line when computing L (FIR) because the emis- 
sion in this line may be contaminated by cloud emission. The es- 
timated values of cooling in all the observed species ([0 1], [C 11], 
OH, H2O, CO) and the total FIR line cooling are summarised in 
Tab.|6] 

The line cooling at different evolutionary stages is shown in 
Fig.|9]by means of histograms of L [Oi], L mo /, and L (FIR)/L/, / 



for Class 0, I, and II sources. The luminosity values reported 
in the figure are for: the Class sources ob served with I SO by 
lGianniniet"ai] d2001l) and with Herschel by INisini etal] d201fl 
(17 Class sources); the Class I (or unresolved Class II + I) 
sources analysed in this paper (T Tau, DG Tau B, and FS Tau 
A+B) , complemented by ISO observations of Class I sources 
bv INisini et al.1 d2002h. a nd PACS observations of HH 46 by 
van Kempen et al. (2010) (14 Class I sources); the Class II 
sources analysed in this work (DG Tau A and RW Aur) (2 Class 
II sources). Note that only 16 out of 17 Class sources and 11 
out of 14 Class I sources are reported in the histogram of the 
molecular cooling ( middle pane l of Fig . [9 V The other sources 
in the sample by Giann inTet alJ (|2001) and INisini et alJ ([2002) 
do not show molecular line emission, probably due to the lim- 
ited sensitivity of ISO. 



Despite the small statistical samples, in particular for Class 
II sources, the histograms in Fig. [9] and the range of values re- 
ported in Tab.|6]indicate that the total FIR cooling decreases with 
the source evolutionary stage going from values of ~10" 3 - 3 L Q 
for Class sources to values of ~ 5 10" 4 - 10" 1 L Q and ~ 9 10" 4 
- 2 10" 3 L for Class I and Class II sources. In particular, L [Oi] 
is of ~ 10" 3 - 4 10" 1 L in Class sources, of ~ 3 10" 4 - 10" 1 
L in Class I, and ~ 10" 4 - 10" 3 L in Class II sources. The 
cooling in the molecular lines, L mo /, shows a stronger decrease 
going from Class (L OTO / up to 2.4 L ) to Class I (L mo / up to 
0.04 L ), and Class II sources (L mo ; ~ 7 - 8 10" 4 L Q ). This is 
due to progressive clearing of the circumstellar material which 
is accreted or transported away by the observed jets. Finally, the 
bottom panel of Fig.|9]shows that also L (FIR)/L,i, / is decreasing 
going from Class 0, to Class I, and II sources indicating that the 
outflow efficiency in radiating away the source bolometric lumi- 
nosity is decreasing with its evolutionary state. Similar results 
are obtained by Karska et al. {in preparation) for a sample of 
Class and I sources observed with Herschel/PACS. 
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Table 6. Cooling in all the detected FIR lines. 





Class 


U |W I J 


L [C n] 


L (OH) 


1 CHiOl 
l t,n/>w; 




L (CO) 




L (FIR)" 








(L ) 


(L ) 


(L Q ) 


(L ) 




(Lq) 




(L ) 




TTau 


11 + I 


1.2 


4.6 10~ 4 


1.3 10" J 


1.5 10~ 3 - 1.6 10" 


-i 


1.2 10~ 3 -9.7 10" 


i 


1.7 10^-4.0 10" 


-i 


DG Tau A 


II 


1.1 10~ 3 


1.8 10~ 4 


4.9 10~ 5 


2.1 10 -5 -4.2 10 


-4 


5.9 10~ 5 - 2.5 10" 


-4 


1.5 10" 3 -2.0 10" 


-3 


DG Tau B 


I 


4.6 10" 4 


1.5 1(T 5 


9.0 10~ 6 


2.7 10~ 6 - 7.0 10 


-5 


8.0 10~ 6 - 4.2 10" 


-5 


4.9 10~ 4 - 6.0 10" 


-4 


FS Tau A+B 


11 + I 


3.3 1(T 4 


2.6 10~ 5 


2.2 10 -5 


3.9 10~ 5 - 9.7 10" 


-4 


2.2 10~ 5 - 5.8 10" 


-4 


4.4 10~ 4 - 1.9 10" 


-3 


RW Aur 


II 


1.3 10~ 4 


0.0 10~ 5 


1.5 10~ 5 


3.1 10 5 - 5.1 10 


-4 


7.4 10~ 6 - 3.0 10" 


4 


1.8 10- 4 -9.5 10" 


-4 


Class * 




1 10 3 - 4 10 1 


3 10~ 4 - 2 10-' 


0-3 10-' 


0-1.2 




- 9 10- 1 




1 10~ 3 - 2.8 




Class I c 




3 1(T 4 - 1 10-' 


7 10~ 4 - 1 10 1 


- 7 10- 3 


- 2 10~ 2 




- 4 10~ 2 




5 10~ 4 - 1.4 10- 




Class II d 




1 1(T 4 - 1 10" 3 


- 2 10~ 4 


1 10~ 5 - 5 10" 5 


4 10~ 4 - 5 10~ 4 




2 10- 4 - 3 10~ 4 




9 10~ 4 - 2 10~ 3 





a L (FIR) = L [Oi] + L (OH) + L (H 2 0) + L (CO) 

* The luminosity values for Class sources are from Gia nnini et ID d2001l) and lNisini etafl d2010T) (17 s ources) 

c The luminosity v alues for Class I sources are from this work (T Tau, DG Tau B, and FS Tau A+B). Ivan Kempen et all d2010h (HH 46), and 
iNisini et "aIll l2002T) (14 sources) 

d The luminosity values for Class II sources are from this work (DG Tau A, and RW Aur) (2 sources) 



Table 7. Mass loss rates derived from the luminosity of the [O i] 63 fim line are compared with Mj et values from optical forbidden 
lines and M acc estimates obtained from UV veiling and/or optical/NIR HI lines. 



Source 


M jet ([Oi] 63) 


M jet (opt) " 




Mjet/Macc 




(M yr-') 


(M yr-') 


(M G yr-') 




TTau 


1.2 10- 6 


1-7 10-' + ? 


0.3-1.5 10-' + ? 




DG Tau A 


1.1 10~ 7 


3 10- 8 - 3 10- 7 


0.5-2 10- 6 


0.05 - 0.2 


DG Tau B 


4.4 10~ 8 


7 10-" 


2.2 10~ 7 


0.2 


FS Tau A+B 


3.1 10~ 8 


? + <2.5 10-' 


2-3 10^ 7 


0.17 


RW Aur 


1.2 10~ 8 


2 10-' - 2-3 10~ 8 


0.034-1.6 10^ 6 


0.008 - 0.35 


Class c 


1 10-' -410- s 








Class I d 


3 10~ 8 - 1 10~ 5 


io- 8 - 10~ 6 


10~ 6 - io- 5 




Class II " 


1 10~ 8 - i io- 7 


10-" - 3 10~ 7 


10" 10 -4 10- 6 





" M ie t (opt) estimated: from optical [Oi], [Sn] line luminosity bv lWhite & Hillenbra nd (2004) (T Tau N, FS Tau B, RW Aur). lHartigan et al.l 
d!995h (DG Tau A , RW Aurl. lHerbst et alldl997l) (T Tau); and from jet density, velocity and radius estimates bv lCoffev et"aT] j2008l) (DG Tau A). 
iPodio et al.l d20TTh (D G Tau B~ andlMelnikoy et al.l f2009l) (RW Aur). 

h estimates by lHartigan et all dl995l): lHartigan & Kenvonl d2003h : iGullbring et"all i2000h : IWhite & Ghe3 d200lh : IWhite & Hillenbrand 

( !2004l) : ICalvetetal.| j2004h : lBecketalJ ( l201(]|) , , . 

c Mjet ([O i] 63) for Class sources are from lGiannini et all d200 lh (17 sources) 

d M iel ([O i] 63) for Class I sources are from this work (T Tau , DG Tau B, and FS Tau A+B),lvan Kempen et al.l ( 120101) (H H 46), andlNisini et afl 
(2002) (14 sources), M jet (opt) are from lHartigan et alj jl994) ; lBacciotti & Eisloffell ( fl999l) : lPodio et al.l d2006T) . M acc fro mlHartigan et al.Ul994l) 
" Mjet ([Oi] 63) for Class II sources are from this work (DG Tau A, and RW Aur) (2 sources), Mj ct (opt) are from Hart igan et aT] {l995); 
ICoffev et alj ( 120081) , M acc from lHartigan et al.1 jl995l) : lGullbring et al.1 jl998l) 



4.5. Mass loss rates 

The mass loss rate, Mj e t, is estimated f rom the [O i ] 63 p lu- 
minosity by using the relationship by iHollenbachl d 1985b . from 
which it is shown that if the ejected material is moving fast 
enough to produce a dissociative J-shock, then [O i] emission 
will be the dominant coolant in the postshock gas for tempera- 
tures of 100-5000 K. Thus, the [Oi] luminosity is a direct tracer 
of the mass flow into the shock, and hence of the mass loss rate, 



jet- 



MjettMoyr- 1 ) = 10- 4 L [Oi]63^m(L ) 



(1) 



This is a si mpler method to estimate M\ et than the use of or. 



tical lines (e.g.. lHartigan et al.ll995l:lBacciotti & Eisloffell 1 99 1 



IPodio et"ai] 2006) because it does not depend on estimates of the 
visual extinction, the inclination of the system, or the geome- 
try of the outflow. On the other hand, the derived estimates are 
based on the assumption that all the ejected material is J-shocked 



and that all the observed [Oi] emission is produced by shocks. 
Thus, if part of the observed [O i] 63 /mi emission arises from 
a photodissociation region due to, e.g., the UV-illuminated out- 
flow cavities, and/or from the disk we may overestimate the mass 
loss rate. 



The mass loss rates derived from the [O i] 63 //m line are 
summarised in Table [7] and compared with mass loss and mass 
accretion rates derived from optical lines and U V veiling (e.g., 
lHartigan et al J 1 995t IGullbring et al J 1 998l 120001) . The mass loss 
rates derived from the [Oi] 63 fim line are also compared 
with the values estimated previously f or Class and I sources 
dGiannini et alJ200lLlNisini et alJ2002h (see top panel of Fig.©. 
In agreement with the estimated FIR cooling, the mass outflow 
rates derived from the [O i] 63 //m line decrease as the driving 
source evolves from values of ~ 1 10~ 7 - 4 10~ 5 Moyr 1 for 
Class 0, to ~ 3 10~ 8 - 1 10~ 5 M @ yr l for Class I, and down to 
~ 10~ 8 - 10~ 7 M yr~' for Class II sources. 
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L(FIR)/L 



-2 

logx 



Fig. 9. Histograms of L [O i] (upper panel), L mo / (middle panel), 
and L (FIR)/L/, / (bottom panel) for Class 0, I, and II sources 
(solid, dotted, and dashed lines). 



The comparison of Mj et ([Oi] 63 pm) with Mj et andM acc val- 
ues derived from observations at optical and UV wavelengths is 
difficult because the latter show discrepancies up to one order of 
magnitude depending on the adopted method. This is mainly due 
to the fact that Mj et and M acc values derived from optical lines 
and UV veiling are highly dependent on the estimates of the 
jet radius and visual extinction (see, e.g., Hartigan et alj fl995; 
iBacciotti & Eisloff elll 1 999t iPodio et al.ll2006f for a discussion of 
diffe rent methods to d erive Mj et from optical forbidden lines and 
Gullbring et al. 1998 for a discussion of the uncertainties affect- 
ing M acc estimates). Thus, for Mj et (opt) and M acc estimates we 
report in Table [7] a range covering all the different values found 
in the literature. However, both for Class II and Class I sources 
the mass loss rate derived from the [O i] 63 pm line luminos- 
ity is larger than or comparable to the maximum value obtained 
from optical and NIR forbidden lines. This suggests that the 
mass loss rate can actually be larger than previously thought and 
the ejection rate can be up to a few percent of the accretion rate 
(Mj et /M acc up to 0.35). This is even more evident if we consider 
that the lower M^/ Macc. ratios are derived by using the high M acc 
values estimate d by Hartigan et al.l d 19951) . More recent work by 
Gullbring et al.l (I 1998b showed that these v alues can be o veres- 
timated by one to two orders of magnitude . ICabritl (120071) found 
similar high Mj et /M acc when considering the most recent and ac- 
curate M acc and Mj et estimates and showed that high Mj et /M acc 
values may have important implications for jet launching mod- 
els. For example, they show that st ellar winds cann ot produce 
such high mass loss rates while X- (IShu et al.lll994l) and Disk- 
dFerreira et al.i r2006) wind models may provide mass ejection to 
mass accretion ratios up to 0. 1-0.25. 



Another fundamental issue is to understand how the mass 
ejection to mass accretion ratio evolves with the source evolu- 
tionary state. The values derived for the Class I and II sources 
in our sample seems to suggest that the Mj et /M acc ratio remains 
constant. However, it is not possible to draw a firm conclusion 
given the large uncertainties affecting the mass accretion rate es- 
timates and the small size of the considered sample. 



5. Conclusions 

In this paper we have analysed Herschel/PACS integral-field 
spectroscopic observations of Class I and II sources in Taurus 
which are known to drive bright optical jets. Thanks to the 
Herschel sensitivity (100-1000 times larger than ISO) we are 
able to detect the FIR counterpart of optical jets from the se- 
lected Class I and II sources for the first time. An exception is 
T Tau which is a bright multiple system unresolved with PACS, 
consisting of a Class II source and a Class I binary system, and 
assoc iated with at least two jets, which has been observed with 
ISO (Spinoglio et al. 2000). We investigate the origin of the de- 
tected atomic and molecular lines by carefully evaluating the 
spatial distribution of the emission on the PACS detector and 
by comparing line fluxes and ratios with predictions from disk 
and shock models. The results of our analysis are summarised 
below: 

- the emission in the atomic [Oi] and [C n] lines is extended 
and spatially correlated with the optical jet emission. In two 
cases (DG Tau B and RW Aur) we also detect a consistent 
offset in velocity in all the spaxels where the [Oi] 63 pm 
line is detected which indicates a gas velocity in agree- 
ment with the values measured for the associated optical jets. 

- the emission in the molecular HjO, CO, and OH lines is 
spatially and spectrally unresolved. However, by using the 
DENT grid of models we show that for typical low mass 
YSO and T Tauri star parameters the irradiated disk surface 
is unlikely to produce the observed large H2O, CO fluxes 
(up to 10~ 16 and 10~ 15 WirT 2 , respectively) even when 
the source is associated with a strong X-ray field. Slow C- 
and J- shocks (V shock ^40 kms -1 and V shock ^30 kms -1 , 
respectively), on the other hand, can reproduce the observed 
line fluxes for an emitting area of diameter of a few tens to a 
few hundreds of AU. Thus, a shock origin is favoured. 

- high- J CO lines (up to CO J— 36-35) and HjO lines from 
high excitation levels (up to E up ~107 K) are detected 
simila rly to what was o bserve d by Ivan Kem pen et al. 
(1201 Oh and iHerczeg et all (120121) for Class and Class 
I outflow sources (NGC 1333 IRAS 4B and HH 46/47, 
respectively). This suggests that lines from high excita- 
tion levels can be shock excited if the density is high enough. 

- the extended atomic emission may be produced by fast 
J-shocks. Shocks with veloc ities higher than 30 km s~' 
with a radiative precursor dHollenbach & McK ee 1989) 
strongly dissociate and ionize the gas giving rise to high 
[0 1] and [C 11] line fluxes, in agreement with the observed 
line ratios ([Oi] 63/145-15-30, [Cn]/ [Oi]<0.17). Excess 
[C n] emission may be due to UV-heated gas in the outflow 
cavity walls. 

- molecular emission may originate instead in 
slow C- or J- shocks, which preserve molecules 
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(Fl ower & Pineau Des Foretsl l2010h . High pre-shock 
densities are required to populate the high excitation H2O 
levels and reproduce the observed line ratios (i.e. J-shocks 
with n~10 4 -10 5 cirT 3 or C-shocks with n~10 6 -10 7 cirT 3 ). 
We cannot exclude, however, that the disk and the warm gas 
in the outflow cavity walls are contributing to the observed 
emission. 

- the cooling in the FIR lines is decreasing as the source 
evolves: for the Class II sources in our sample the cooling 
is from one to four orders of magnitude lower than for 
Class I and sources (L [Oi] -10 4 -10 3 L , L H 2 and 
L CO ~10~ 4 L Q ). The molecular cooling is decreasing more 
abruptly as the source evolves indicating that for Class 
sources the main coolants are water and CO, while in Class 
I and II [0 1] becomes an important coolant. 

- the mass loss rate for the Class II sources in our sample 
is up to three orders of magnitude lower than for Class 
and I sources, i.e. M jet ([0 1] 63 ^m) ~ 1CT 8 - 1(T 7 M G yr" 1 . 

- the mass loss rates inferred from the [Oi] 63 um line are 
larger than or comparable to values obtained from optical 
and NIR forbidden lines, implying higher mass ejection to 
mass accretion ratios, up to 0.35. This may have important 
implication for jet launching models. 

The above summary places the FIR emission from Class 
II and I jet sources within an evolutionary picture. The Taurus 
optical-jet-sources studied in this work show FIR atomic and 
molecular emission similar to that previously observed with ISO 
for Class and Class I sources, including a highly excited molec- 
ular component. However, the emission associated with Class II 
sources is fainter and more compact (in particular the molecu- 
lar component), and the FIR line cooling and mass loss rates are 
one to three orders of magnitude lower than those estimated for 
Class and I sources. 
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Table A.l. Target, observational mode (LineSpec or 
RangeSpec), operational day (OD), and identifier (OBSID) of 
the analysed HerschelfPACS observations. 



Target 


Obs. Mode 


OD 


OBSID 


TTau 


LineSpec 


272 


1342190353 




RangeSpec 


272 


1342190352 


DG Tau A 


LineSpec 


273 


1342190382 




RangeSpec 


273 


1342190383 


DG Tau B 


LineSpec 


316 


1342192798 




RangeSpec 


678 


1342216652 


FS Tau A+B 


LineSpec 


316 


1342192791 




RangeSpec 


641 


1342214358 


RW Aur 


LineSpec 


290 


1342191359 




RangeSpec 


290 


1342191358 
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Appendix A: Herschel/PACS observation identifiers 



Table lA.ll list the targets, observational modes, operational days 
(ODs), and identifiers (OBSIDs) of the Herschel/PACS observa- 
tions analysed in the paper. 

Appendix B: Line extended emission 
characterisation 

In this Appendix we describe the procedure which is applied 
to ascertain whether the line emission is extended and/or offset 
with respect to the continuum emission. If both line and contin- 
uum emission originate from the same region, supposedly the 
star-disk system which is unresolved with PACS, any emission 
detected out of the central spaxel is due to the spectroscopic PSF 
and/or to the fact that the source is not centred on the central 
spaxel. The line and the continuum PSF have the same shape and 
are centred at the same position, hence the line-to-continuum 
ratio is constant across the PACS field of view. If, on the con- 
trary, the line emitting region is more extended and/or offset with 
respect to the continuum emitting region the line and contin- 
uum distribution across the PACS field of view are different and 
the line-to-continuum ratio should vary in the different spaxels. 
However, when most of the line emission is emitted close to the 
source it can be difficult to detect extended and/or offset emis- 
sion given the low spatial resolution and sampling offered by 
PACS. 

To check for the presence of extended line emission we sub- 
tract from the line + continuum image the "on-source" emission 
and search for "residual" local emission above the confidence 
level. 

An image of the line and continuum emission is con- 
structed integrating the PACS data cube over the wavelength 
range covered by the considered line. For example to obtain 
line+continuum image for the [Oi] 63 yum line we integrated 
on the spectral range from 63.141 fim to 63.242 fim, for all 
the sources in our sample. Integrated fluxes, F/ ine+fo „,, and er- 
rors, &Fime+cont> m units of W irT 2 ', have been computed for each 



spaxel summing the flux density on each spectral element, fi, in 
Jy, over the n spectral elements in the defined wavelength range 
as follows: 



line+cont 



= dvY j f i 



AF, 



line+cont 



— dv - 



1 !=1 



(B.l) 



(B.2) 



where dv is the average spectral element size, in Hz. The 
continuum image, F conl , with the associated error, AF co ,„, in Jy, 
are estimated in each spaxel by computing a weighted average 
of the flux density over a region of nl spectral elements adjacent 
to the detected line: 



Wi = 



1 

A/f 



AF r , 




(B.3) 



(B.4) 



(B.5) 



Then we subtract from the line + continuum image, F;,„ e+fo „,, 
the "on-source" emission, i.e. the line + continuum flux in the 
spaxel showing the brightest continuum, which is scaled accord- 
ing to the continuum level in each spaxel. If j is the spaxel where 
the continuum emission is maximum, we obtain the image of the 
residual flux, F res id m ,i, and the associated error, A¥ res id ua i, as: 
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FcontJ 



me+cont,j 



ARj = Rj 
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The confidence level at which residual emission is detected 
in each spaxels is: 

F residual /T ~ t , 

^A-F-d-i (R10) 

Fig. |2] (right panels) shows the displacement of the [Oi] 63 
fim residual emission, F resu j m i, with respect to the continuum 
emission, F cont , and the optical jet direction. We clearly see an 
offset between the continuum emission (dotted contours) and the 
[O i] 63 fim line residual emission (solid contours) which is dis- 
placed along the optical jet PA (blue/red dashed lines). The spa- 
tial correlation between the [O i] 63 fim line residual emission 
and the optical jet is evident also for the sources for which the 
line and continuum emission peak on the same spaxel before the 
subtraction of "on-source" emission (e.g., T Tau and RW Aur; 
see, for comparison the maps in the left panels of Fig. 0. 
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